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PREFACE 


The purpose of this research was to examine quantitatively the 
possible effect of the slope of the core of a cyclonic system on the 
pressure tendency at the center and the further relationship, to that 
same tendency, of various other synoptic values available to the fore¬ 
caster. 

This work was conducted at the U. S. Naval Postgraduate School, 
Monterey, California, during the period November 1948 to May 1949. 

I wish to express my appreciation to Professor William D* Duthie 
and Associate Professor A. B. Mewborn for their advice and guidance 
throughout. 
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INTRODUCTION 


The variation of pressure with time in any meteorological situation 
has long been recognized as one of the basic foundations upon which 
prognosis and forecast must rely* As a result, a rather extensive amount 
of meteorological literature has been devoted to the tendency equation in 
its various forms, and possible kinematic applications thereof* Prom a 
qualitative standpoint, the theory has been well developed* Quantitative 
studies have, however, progressed more slowly under a number of widely 
recognized handicaps. Although the principal work has been done on the 
local pressure variation, e.g. the recent work of Godson (6), certain of 
the rules developed, such as that of Petterssen (11), are concerned with 
the deepening or filling of the center of the system itself. These re¬ 
sults require either moderate to extended computations, data not ordinar¬ 
ily available under routine conditions, or are applicable, with a 
sufficient degree of accuracy, for only relatively short range forecasts* 
The practicing meteorologist, under pressure of a fixed forecast deadline 
and with a forecast period of at least twenty-four hours, finds that he 
must still rely on experience and theoretical judgement in determining 
the probable pressure variation to be expected of a given system. 

This study was undertaken to ascertain whether some method of 
quantitative assistance could not be uncovered through statistical in¬ 
quiry. In order that any favorable results might be of use to the prac¬ 
ticing forecaster, the following conditions were established: 

(a) Availability of Data - That such data as might be necessary 
for the obtaining of the desired results should be included in that 
presented by the several maps customarily analyzed in the average 
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meteorological office* Further, that it should not be limited by being 
among that data received from only selected stations, i .e. , radiosonde 
or pibal reports. 

(b) Ease of Application - That the results be adaptable to repre¬ 
sentation in such a form that extended computations would not be 
necessary. 

(c) Extended Period of Application - That any results be valid for 
a period of at least twenty-four hours since that is the usual prognostic 
interval. 

With these limitations in mind, and in view of the fact that the 
future intensity of a particular cyclone is ordinarily one of the 
major items of interest to the "customer" of the meteorologist, the 
investigation was confined to the future deepening or filling of the 
centers of low pressure systems. 

For certain theoretical reasons, the original study was directed 
towards a determination of the possible correlation between such deep¬ 
ening or filling and the slope of the core of the cyclonic center as 
expressed by the horizontal distance between the surface center and 
the corresponding centers at the 700mb and 85Qmb levels. The assump¬ 
tion was made that this core was a straight line, at least from the 
surface to the 700mb level. Although suoh an assumption was believed 
to be incorrect, it was felt that the error introduced by variation 
from that condition might prove to be sufficiently small as to be 
negligible. Such as shown to be the case. The investigation was 
later expanded so as to inquire into the possible individual or 
further correlation of the following additional factors* 
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(a) The pressure at the center of the system* 

(b) The speed of the system* 

(c) The E-W component of the speed* 

(d) The N-S component of the speed* 

(e) The horizontal distance between the actual position of the 
center at the 850mb level and the position indicated by the assumption 
of a straight core from the surface to 700mb* 

(f) The deepening or filling of the system during the previous 
twenty-four hours* 

It soon became apparent to the author that a certain mutual corre¬ 
lation existed between the prospective deepening or filling of the 
system, the originally established criterion, and the existing pressure 
of the cyclone’s center* At this time (b) - (f), above, were re-examined 
to determine their possible connection with the discrepancies between 
the calculated pressure change, during the ensuing twenty-four hour 
period, and that observed to occur* 

As it developed, the investigation, based originally on the work 
of Bjerknes and Holmboe (2), divided itself into four parts as follows* 

(1) To ascertain whether there existed any particular distance 
and hence any critical slope angle of the core, which would serve as 
a line of demarcation between future negative and positive pressure 
changes* 

(2) To determine whether there existed any easily expressed re¬ 
lationship between the slope, expressed as a distance, and the magnitude 
of the pressure change* 

(3) If (2) was affirmatively concluded, whether any of the other 
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variables would satisfactorily explain the remaining discrepancies* 

(4) And finally, if (3) yielded positive results, whether any 
third variable could account for the residual error* 

Since this was a statistical investigation, the resulting re¬ 
lationships are presented in terns of simple and multiple correlation 
coefficients, two, three, and four-variate regression equations and 
standard errors of estimate* 

Data was assembled from maps covering a period of six consecu¬ 
tive winter months, October to March inclusive. Map time for the 
surface analysis was 0030Z, for the upper levels 0330Z, and for the 
limits of the twenty-four hour period 0030Z; the reasons for such 
choice of times will be discussed later. Sufficient cases were in¬ 
cluded so as to render the results, if not conclusive, at least 
reasonably reliable on the basis of standard sampling theory* In 
view of the intention to limit all computations and data to within 
that reasonably to be expected of a practicing aerologist, a maximum 
of three-figure accuracy was used — two-figure where the probable 
accuracy of analysis and visual subdivision of measurements would 
so indicate. 

With the above considerations in view, the following results 
of the investigation are advanced; 

(1) That, within rather small limits of error, there does 
exist a critical value of the horizontal distance between the centers 
of the cyclone on -the surface and at the 70Qmb level with respect 
to the future difference of pressure at the center of the system* 

This value for Di is 2,83°Iat. or 170 miles. For Dj_ less than that 
value, the system may be expected to have filled at the end of the 
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next twenty-four hours for D]_ greater them 170 miles to have deepened. 

(2) That the correlation coefficient between Di and dp is large 
and negative, being -.70 and that the multiple correlation coefficient 
of p p based on Dj and the original pressure at the center of the system 
is large and positive, .73 * These correlations, especially the latter, 
indicate a pronounced relationship among the three variates and a reason¬ 
able degree of reliability might be expected of the regression coefficients 
computed therefrom. The corresponding standard errors of estimate are 
6.34 mb and 5*2 mb and these are rather large in view of the mean value of 
dp, which is -4 mb. However, since dp varied from a+12mb to -25 mb, a 
range of 37 mb, the latter standard error of estimate defines a limit of 
25-30/S of the range. It might be pointed out that a more proper com¬ 
parison may be made by the use of actual pressure rather thandp. In 
this case, since the extremes of surface pressure may be considered to 

be 850 mb and 1050 mb, a dispersion of 200 mb, either Sp represents an 
error of only about 3«0^. 

(3) That, through the use of the computed regression equations or 
the corresponding graph, the analyst may reasonably expect to forecast 
the central pressure of a cyclone twenty-four hours in advance with the 
probability being .50 that his predicted pressure will be within 5 mb 
of that observed at the end of the period. It is believed that, within 
this not unreasonable limitation, such a result is of some assistance in 
quantitative forecasting. 
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CHAPTER I 


EARLIER INVESTIGATIONS OF PRESSURE TENDENCY 

1* Theoretical derivations. 

The first, and what still might be called the classical, devel¬ 
opment of the tendency equation was carried out by Margules (?) in 
1?04* The later treatment by Bjerknes (1) and Holmboe (7) differ only 
in notation. Beginning with the hydrostatic equation and advancing 
through familiar steps, we arrive at the equally familiars 



which, for the surface pressure tendency, reduces to 
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on the recognition of the fact that there can be no vertical motion 
through the level of the surface of the earth. 

On the same basis, but through a different evaluation of the 
hydrostatic equation, Houghton and Austin arrived at the following ex¬ 
pression for the surface pressure tendency* 
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in which, they have shewn, the two terms on the right side of the 
equation are of the same order of magnitude. 


An adaption by Raethjen (12) of an original derivation by 
Panofsky (10) provides us with a third equation. This development is 


based on the series of identities: 


(t) 
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and results in the following expression for the surface pressure 
tendency: 
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where/ 1 is the dry adiabatic lapse rate,o** the actual lapse rate, and 
the e corresponds to an expression involving 2^* and represents the 
results of non-adiabatic temperature changes. 

Recent work by Godson (6) has produced still another result which 
he chooses to call the "isobario tendency equation"* Here, from the 
mathematical sequence of : 


y> _ /x \ _ z' 

*2>£- ^ C 2>z? ' S? ^ / 

- ^2T je ^7 ~ 2>/> 

~r 7>i P 2>* 2>£r 
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he arrives at an expression for the surface pressure tendency as 
follows: 
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2* Quantitative approaches* 

There have been at least three proposals for the reduction of 
pressure tendency to a calculable status: 

(a) The technique of Rossby as applied to the "classical” 
equation after the inclusion of certain limiting assumptions by 
Bjerknes and Holmboe (2). 

(b) Godson's (6) evaluations of his "isobaric" tendency 
equation* 

(c) Certain rules developed by Petterssen (11) from the gen¬ 
eral mathematical expression for the difference between the "local” 
sind ’’individual” changes of any variable; in this case, pressure* 

All of these lend themselves to practical application, to a vary¬ 
ing degree. However, each fails, in at least one respect, to meet the 
requirements, as established for the purpose of this study, of a practical 
forecasting tool. A brief discussion of the details of each technique 
follows - and from this it will be seen that Petterssen's miles come 
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the closest to fulfilling the desirable qualifications referred to in 
the Introduction. 

The technique referred to in (a) is not strictly applicable to the 
problem under discussion, since it is concerned with the determination 
of the pressure tendency at some upper level, based on the tendency at 
the ground. However, it does serve as an example of early efforts to 
correct the tendency equation into a form susceptible to mechanical com¬ 
putations. In reference to the "classical" tendency equation: 



Bjerknes and Holmboe (2) attempted to show that, with certain 
assumptions, the last two terms on the right side were of opposite sign 
and that their joint contributions could be neglected reducing (8) to the 
one or "advective" term. Rossby then derived from this sin equation for 
quantitative solution of pressure tendency. In so doing he introduced 
further assumptions, among them that of geostrophic wind. Unfortunately, 
both of these developments were later shown to be invalid. However, even 
if valid, Rossby's equation would not have been acceptable under the con¬ 
ditions established for this investigation since it involved the area 
under a wind-shear hodograph from 3,000 - 10,000 feet. This restricted 
its application to stations where such upper air soundings were available. 
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Godson (6) gives an example of the computation of the surface 
tendency based on the observed tendency at 15,000 feet. From 
equation (7), he develops the applicable form to be: 



In evaluating the several terms, he finds it necessary to use 

a) A graph of values of P^, observed at two-hour intervals, 

b) The wind hodograph from the friction layer to 15,000 feet 
corrected for non-geostraphic wind values, 

c) The appropriate tephigram, and 

d) Elsasser’s radiation chart 

with supplementary computations in each case. The necessity- for such 
data and calculations, in addition to its short period of reliability, 
disqualifies this technique under that same set of limitations advanced 
earlier. 

We arrive, then, at Peterssen’s (11) adaption of equation (2) to 
quantitative use by its utilization only at points where the convective 
tem(£-V/j)vanishes. While it is undoubtedly true that, at such points, 
the instantaneous ’’local” and "individual" changes are equal, that is 



the extension of the time interval introduces the obvious fallacy of 
strict adherence to extrapolation* 

It is recognized that no claim was made in either of the first two 
oases other than as supplementary computations leading to rough first- 
approximations in support of dynamic theory, nor on the part of the 
third to other than qualitative interpretation* 
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CHAPTER II 


DETAILS OF THE INVESTIGATION 

1* Organization of data* 

The period covered by the investigation was from October, 1945 
through March, 1946. The charts used were the 0050 surface (sea level) 
and the 0550 upper air (700 mb and 850 mb constant pressure)* These 
charts were obtained from the chart file of the Department of Aerology, 

U. S* Naval Postgraduate School, Monterey, California. In order to com¬ 
pensate for the time difference between the two maps, the position of the 
surface center was advanced one-seventh of its observed travel during the 
following twenty-four hour period. This was done in the interest of in¬ 
vestigative procedure although it is the author’s belief, as will be dis¬ 
cussed, that this was not necessary. 

In selecting the cyclones to be used in the investigation the follow¬ 
ing precautions were observed? 

a) The charts ware carefully examined with the purpose of ensuring 
the correctness of the analysis. 

b) The surface centers were enclosed by at least one isobar and 
were sufficiently well defined as to permit the location of the center with 
reasonable accuracy. 

o) The centers on the upper air charts were, in the majority of 
cases, enclosed by one or more isohypses. However, in a few instances, 
where an abnormally wide spacing of isohypses and one or more reports per¬ 
mitted the reasonable location of a center not encircled by a standard 
isphypse, such a center was used. 
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d) In general, only those areas of the charts were used where 
adequate reports guaranteed the certainty of the analysis in both the 
surface and upper air charts. However, in order to test the extension 
of the theory to two areas of meteorological importance not meeting this 
qualification, a number of cases in the Aleutian Islands-Alaskan Gulf and 
Hewfoundland-Nova Scotia regions were included. 

A total of eighty-four surface cases were tested. Of thesej seventy- 
five had a corresponding center on the 700 mb chart, seventy had a corres¬ 
ponding center on the 850 mb chart, and sixty-one had centers on both 
levels, associated with the particular surface cyclone. The following 
data, where applicable, were recorded or computed in each case, with all 
distances being expressed in degrees of latitude: 

a) date-time group (for identification) 

b) surface center - latitude 

c) surface center - longitude 

d) surface center - central pressure 

e) surface center - V x 

(Since the distance was expected in the majority of cases to be to the west, 
this direction was taken as positive for all latitudinal velocities and dis¬ 
tances — despite its contradiction of the accepted custom). 

f) surface center - Vy 

g) surface center - corrected latitude 

h) surface center - corrected longitude 

i) 850 mb level center - latitude 

j) 850 mb level center - longitude 

k) 850 mb level center - V x 
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l) 850 mb level center - Vy 

m) 700 mb level center - latitude 

n) 700 mb level center - longitude 

o) 700 mb level center - V x 

p) 700 mb level center - Vy 

q) 850 mb to surface - difference of latitude 

r) 850 mb to surface - difference of longitude 

s) 700 mb to surface - difference of latitudes 

t) 700 mb to surface - difference of longitude 

u) horizontal distance - 850 mb to surface 

v) horizontal distance - 700 mb to surface 

w) pressure change - next twenty-four hours 

From among these, or items computed from them, were obtained the variates 
to be used in the several correlations carried out# 

2, Tetrachoric distributions. 

On the basis of Bjerknes (l) original work on waves in the westerlies, 
a discussion of which may be found in any advanced text, it is reasonable 
to predict a probable marked degree of relationship between the horizontal 
distance of the surface center of a cyclone from the center at any upper 
level; i.e., the slope of the core of the system, and the magnitude and 
sign of the value for Up given by the standard tendency equation# Quali¬ 
tatively, the reasoning may be expressed in this manner. 

a) The divergence term depends on the distribution of divergence 
and convergence above the point under consideration, the line of demar¬ 
cation being the level of the last closed isobar. This level, for a given 
system will have a close relation to the horizontal distance between the 


(14) 



surface and upper centers. From the nature of the dependency, it might 
be expected that a decrease in the magnitude of the distance would re¬ 
sult in an increase in the algebraic value of the pressure tendency. 

b) To a great extent, the advective term will depend on the 
thermal asymmetry or symmetry of the cyclone. Basically, such asymmetry 
will be interconnected with the same distance as above. It might further 
be expected that the nature of the interaction would be the same as in (a). 

On the basis of (a) and (b), the reasonable assumption was made that 
the total pressure tendency at its center would vary directly with the 
slope of the core of a cyclone, or inversely with the horizontal distance 
between the centers of any two layers, and that there might be some critical 
value of the latter at which the corresponding A p would change from negative 
to positive. The layer over tvhich the horizontal distance would be measured 
was chosen as the surface to the 700 mb level, since this was the thickest 
layer for which charts are usually drawn at each boundry and through which 
the center usually maintains its identity. The distance so measured was 
considered as The corresponding values ofAp and were plotted, 

resulting in Fig. 1. 

Since the initial interest was in the possibility of the existence of 
a critical value for Dj, this was determined, by inspection, to be 2.8j° lat. 
(170 mi.) and produced the tetrachoric distribution shown: 


^P,<e.83 a I 0,>^8 3 


4p C+) 
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This distribution was evaluated in two ways. First, it might be 
considered that quadrents II and IV include the theoretically correct 
cases in that they support the initial assumption. On this basis, sixty- 
seven of the seventy-five cases, or 89 » 5 ^> successfully agree with the 
prognostic criterion. It might be pointed out that, if the best weighted 
estimate of quadrent distribution were used} i.e., in the proportions in¬ 
dicated by the end totals of columns and rows, the "correct” percentage 
would be 65 »5^. If a pure chance or equal quadrant distribution were 
assumed, a percentage of 50.0?£ would result. Secondly, the improbability 
of such a distribution occurring by chance was expressed by the Skill Score 

j|( 

Method of Brier (3). The resulting Score was 73»8/». If the above weighted 
and chance distributions were substituted the Skill Scores would be Q,2% 
and 0.0/ respectively, which even more markedly indicates the improvement 
of this distribution over either random or intelligent guesses. 

Fig. 2 is a plot made of corresponding values of dp and D 2 , where 
D 2 represented the horizontal distance between the surface and 850 mb level 
centers. Here no critical distance of demarcation was apparent and no 
tetachoric distribution was possible. Possible explanations for this are 
that the layer is too thin to be indicative of the desired slope or that 
this level might still experience undue influence by surface friotion 
and/or friction layer turbulence. 

3. Two-variate correlations. 

Since the relationship between and dp presented at least a tetra- 
choric distribution, it was reasonable to inquire into the possibility of 

*See Appendix for a description of the method of computation. 
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a quantitative statistical correlation between the two variates. The 
scatter diagram resulting from the plotting of two variables in con¬ 
junction, Fig. 1, indicated a favorable prospect for such a possibility. 
Actual calculations produced the values for the following simple sta¬ 
tistical functions: 

(a) simple correlation coefficient,/*, of the deepening or filling 
to the displacement between the two centers 

(b) the arithmetic means ofap and 

(c) the unit of standard deviation (variance from the mean), <T , and 

(d) standard error of estimate, S. 

(Xi s Apj Xp s Dfc) 

\ a -4.85 (mb) X2 a 4.4? (°lat.) 

*1 = 9.19 = 2 *55 

/*^ 2 —*70 S|2 — 6.52 * 

The regression equation of Xi on Xp as calculated from this data was: 

(12) X x a 2.76 X2 +7.44 

and this regression line is plotted in Fig. 1. As a matter of investigative 
curiosity and, principally, as a preliminary to later multiple variate cor¬ 
relations, a number of further relationships were computed; namely, 

X| -4p - on 

Xt - I >2 - the distance between the system centers at 
y *" surface and 8 50 mb level 

.. Xtj - P - existing central surface pressure 

X0 - V - velocity of storm along its course for succeeding 
twenty-four hour period 

* Details of computations may be found in Appendix 
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X 2 - Di on 
x 5 - P 
Xs - V 
Xg - V - on 
x 5 - P 

The simple statistical functions, as defined above, are listed below 


a) 

n 

- 6.26 

(mb) 


*1- 

8.5? 



^15 * 

-.2? 


*) 

X 2 a 

-5.1? 

(mb) 


^1 = 

?.04 



/*15 a 

-.55 


c) 

% * 

-4.70 

(mb) 


(Ti = 

?.54 



'•a = 

.55 


a) 

x 2 ■ 

4.45 

(°lat.) 


<5~2 a 

2.54 



r 2$ = 

.55 


e) 

X 2 a 

4.25 

(°lat.) 


^2 - 

2.52 



^28 * 

-.55 



Xj > 2.86 (oiat.) 


II 

2.14 



s 15 ■ 

: 8.25 



Xj a 

995*82 

(mb) 


<r 5 « 

11.57 



S15 , 

» 8.50 

♦ 


lc 8 = 

?.45 

(°lat./24 

hr) 

*8 = 

4.25 



s l8 ■ 

i 8.78 



*5 * 

??5.65 

(mb) 


*5* 

11.51 



s 25 * 

. 2.20 



X8 a 

9.55 

(°lat./24 

hr) 

<^8 = 

4.15 



s 28 * 

1 2.18 
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f) Xg = 9.45 (°lat/24 hr) X 5 = 990.25 (mb) 

<Tq = 4.Z5 ^5 » 10.52 

/ '*85 * ”* 12 s 85 * 4 * 2 5 

The scatter diagrams for the above pairs (except (a)) of inter¬ 
related variates are shown on the succeeding pages as Figs. 3-7. 

In addition to these computations, a visual inspection was made of 
additional data to ascertain the possibility of further correlations. The 
following relationships were inspected: 

X I - 4p - with 
X — 

4 (1 - D 2 ) - the horizontal distance between the actual 

2 center at the 850 mb level and the estimated 

center based on the assumption of a straight 
line core to the 700 mb level 

X 6 - V x - the y-component of the velocity of the system during 
the ensuing twenty-four hours 

- the deepening or filling during the previous 
' twenty-four hours 

None of these showed evidence of sufficient correlation to warrant calcu¬ 
lation of the statistical functions. 

4. Multiple-variate correlations. 

In order to complete the statistical investigations a number of 
multi-variate correlations were carried out. These, and their associated 
multiple correlation coefficient, were: 

a) X-^ (dp) on X 2 (Di) and X 5 (P); » .93 * 

b) Xj_ on X 2 and Xg (V); r 1#2 g 8 

c) Xj on X^ and Xg; s .46 

d) Xj on X^,, X^ and Xg; r 1 # 2 g 5 * .77 

* Details of computations may be found in Appendix. 


( 21 ) 



d V 







(25) 






(24) 









The basis of these computations was to determine whether any 
combination of independent variates might interact in such a manner as 
to govern the expected values of the dependent variable more closely. 

With the exception of the original proposal of the existence of a re¬ 
lationship between the deepening or filling of the center of a system 
and the slope of its core, no attempt was made to predict, on the basis 
of theoretical considerations, the possible effects of the other variables 
acting either separately or jointly. This explains the rather random 
choice of variates with which the investigation was carried out - the 
choice being limited, in fact, only by the amount of information available 
to the forecaster on the two maps under consideration. In this connection, 
the question might be raised as to the propriety of using the velocity of 
the system during the twenty-four hours following the map time for the sur¬ 
face analysis, as both a variable and as the guide for the advancing of the 
surface system up to map time for the upper level analysis. Obviously, this 
fails to satisfy the original requirement that variates used as an argument 
in any forecasting technique developed be available to the analyst on the 
two maps utilized and not require prognostication themselves. 

In accordance with statistical practice, a series of independent cases 
were tested to verify the results of the investigation. These will be 
discussed in full in the next chapter, it being sufficient to state at this 
point that in these test cases the velocity for the past twelve hours was 
extrapolated ahead for the three hours lapse between the surfaoe and 700 mb 
maps. The difference in the majority of the cases in the value of dp was 
found to be less than 2.5 mb from the value obtained by use of the future 
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movement for advancing the surface system. In those instances where 
the discrepancy was larger than 2.5 mb; the number resulting in in¬ 
creased accuracy was slightly greater than those introducing added 
inaccuracy. The actual figures were as follows: 

Number of random test cases 50 

Number in which discrepancy <2 ,5 mb 31 

Number in which discrepancy>2.5 mb 
but resulted in closer approach 
to observed (dp) 11 

Number in which discrepancy >2*5 mb 

and introduced less accurate result 8 
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CHAPTER III 


SUMMARY OF CONCLUSIONS 

1* Evaluation of results* 

The principal conclusion that may be drawn from this study is the 
sustaining of the assumption that there is a high degree of correlation 
between the slope of the core of a cyclonic system and its future in¬ 
tensity trend — and, further, that some critical value exists at 
which this trend changes from negative to positive. This latter may 
be substantiated from the evaluation of tetrachoric distribution of 
the previous chapter. A percentage accuracy of 89.5$ and a Skill Score 
of 73.8/6 are sufficiently large to establish the credibility of the 
distribution and the existence of a critical distance of 2*83° latitude 
or 170 miles. With a distance of less than this amount, the system 
may, with a great deal of reliability, be predicted to fill and, 
conversely, with Dj equal to or greater than that value, to deepen. 

The computed two-variate correlations effectively served to 
establish more completely the inter-relation between the above two 
variables. The correlation coefficient was large and negative 
(-.70) and was markedly greater than the remaining coefficients, 
whose magnitudes changed from .12 to .35. TThile such a value of 
r« -.70 is customarily accepted as the lower limit of the statisti¬ 
cally significant range, it represents in the field of meteorology 
a relatively important correlation, since the numerous intangible 
variables, acting simultaneously, reduce the threshold of credi¬ 
bility to a decidedly lower value. Because the correlation is nega¬ 
tive, a decrease in the displacement will result in an increase in 
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the future pressure. It may be seen from Fig. 1 and the magnitude 
of the units of standard deviation and the standard errors of estimate 
that considerable dispersion from a linear relation still remains. It 
was in an attempt to reduce this dispersion that the multi-variate 
correlations were undertaken. 

With one exception, these multiple relationships did little to 
improve the situation, their magnitudes being .46, .71, and .77. 

Even the latter two values are not sufficient to warrant the added 
consideration of the second variable. However, the exception, which 
was the mutual consideration of the principal distance (from the 
above section), the existing pressure at the center and the inten¬ 
sity trend of the system, resulted in a multiple correlation coe¬ 
fficient of .93. This value is considerably better than any of the 
other relationships investigated and can be considered statistically 
significant without reservation. From this value of r^.25* ' fche 
following regression equation was computed: 

X X - -3.62 X 2 -.52 (X 5 - 994) 11.00 

where Xj, -£p during next twenty-four hours, in mb© 

-^2-^1 displacement between surface and 700 mb. 
centers in degrees latitude. 

and X_ - P, existing central pressure of system, 

5 in mb. 

This equation is readily expressed in graphical array, as in 
Fig. 8, for ease of utilization. By means of this graph, estimated 
values of Xj were obtained for the original data and such estimated 
results were plotted against the observed pressure changes from 
that data. Fig. 9. It can be seen that the correlation is linear. 
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the slope of the regression line equals l.fl and that 5 Of* of the oases 
lie within 5.2 mb. of that regression line. Admittedly, this is 
large compared with the mean pressure change,/! p~ -4, mb., but (as 
was pointed out in the Introduction) compared to the range of either 
pressure change or the surface pressure (of the order of 200 mb) its 
magnitude is not so alarming. Certainly, the obtaining of an esti¬ 
mate of a central pressure, twenty-four hours in advance, which should 
be within 5 mb. half of the time is of some value, especially to the 
inexperienced forecaster. That this reliability is fairly consistent 
may be seen from the fifty supplementary cases. In Fig. 10, these 
are plotted with the same arguments as above. The regression line 
is the same and BOfo of these test cases lie within 4.9 mb. This 
agrees closely with the original results and serves to emphasize 
the limits set by them. 

2. The level of non-divergence. 

According to the classical theory, the level of non-divergence 
coincides with the level of the last closed isobar over any £ven 
point. ' 'When the mass divergence above that level equals the mass 
convergence below, the pressure tendency at the surface will be zero. 
If, in equation (13), , or^p, is set equal to zero and certain 

assumptions are made concerning the intensity and dimension of a 
cyclone, then it is possible to calculate the corresponding height, 
of the last closed isobar above the center. For this purpose, the 
following three cases were assumed* 

Central pressure - 1002 mb. 993 mb. 975 mb. 

Diam. of surface system - 300 mi. 600 mi. 1000 mi. 

Diam of 700 mb. system - same as surface 
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Although the author has no statistics to support them, it is 
felt that these figures are quite representative of three stages of 
cyclonic development. The corresponding levels of non-divergence 
for these cases were computed to be 12,200 ft., 14,900 ft., and 
13,700 ft., respectively. In this connection it might be pointed 
out that Charney (5) has shown that Rossby's formula for the speed 
of a barotropic wave will hold for baroclinic waves if the mean 
zonal wind is used and that this is about equivalent to the wind 
at 600 mb* He further quotes Holmboe to the effect that this is 
true only at the level of non-divergence, which also must bo at 
approximately 600 mb. The data of both Willett (13) and V* 

Bjerknes (3) on the level of average zonal wind, and hence of 
non-divergence, agree quite closely on the value of 580 mb* or 
4*25 km* (13,800 ft*). The close agreement of the computed 
values to the above serves in the nature of substantiation of 
the investigation. It is believed that the variation might have 
been introduced by the assumption that the diameters of the two 
centers are equal* It is the author’s contention that many cases 
of non-conformity to the miles and results of the previous para¬ 
graph are those in which this assumption does not hold* 

3* The forecasting method* 

The procedure for applying the above results needs little ex¬ 
planation. Entering the table of Fig* 8 with the arguments of cen¬ 
tral pressure and horizontal distance between the 700 mb. and sur¬ 
face centers, one obtains the expected pressure change during the 
next twenty-four hour period. Admittedly, there are a number of 
limitations to application of this study. They may be itemized 
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as follows* 

a) The presence of a center on the 700 nib, level is required# 
This need not be a closed center if there is sufficient evidence, of 
the nature of an inconsistent shallowness of pressure gradient,to in¬ 
dicate reliably the existence of an incipient center, 

b) Every effort should be made to ascertain, from the avail¬ 
able reports, the most probable position of the lowest point of the 
700 mb, surface, since reliance on the geometric center of the lowest 
closed isophypse will often introduce error, 

c) Considerable difficulty will be encountered in the case 
of a single upper center associated with two adjoining cyclones on 
the surface. Unless this occurs in a region in which upper-air 
soundings are relatively numerous, a correct solution is improbable. 

To summarize* extreme care must be taken with the 700 mb, analysis 
to insure sufficient accuracy for the application of this Systran, 

It is believed that this forecasting tool will find its greatest 
application to the prediction of pressure tendency of those severe, 
winter storms whioh, originating in Western Canada, recurve in the 
Chicago-New York area, 

4, Possible further investigation, 

A number of avenues of further study are opened by the results 
obtained above. It might prove interesting, for example, to inquire 
into the possibility of a similar relationship existing in the case 
of anticyclones, 

A great deal of further work could be done in connection with 
paragraph 2 above. Actual statistics on the dimensions of the surface 
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and 700 mb* centers in conjunction with the central surface pressure 
and tendency would allow at least three correlations to be investi¬ 
gated t 

a) The calculated horizontal distance between the two centers 
with respect to the observed distance* 

b) The calculated level of non-divergence with respect to the 
observed pressure tendency. 

c) The discrepancy between this calculated level and 13,800 ft* 
with respect to the difference between the computed pressure change 
and that observed* 

Such lines of inquiry are, of course, in addition to the verifi¬ 
cation by trial of these results and their integration into the theory 
of dynamic meteorology* 

Had time permitted the author would have been interested in incorp¬ 
orating the wind and density at stations to the southeast a distance of 
one-half the radius of the surface center in eaoh of the above cases 
into the correlations as further representations of the advective term 
of the tendency equation* 

In conclusion, it is felt that the principal relationship obtained 
by this inquiry, as presented in Fig, 8, is of a sufficient high degree 
of reliability to permit its use as a quantitative forecasting tool, 
bearing in mind the limitations set forth in the preceding section* 
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APPENDIX 


1* Statistical computations# 

All correlation calculations in the paper were based on stan¬ 
dard statistical principles as defined by Kenney (8)# Basically, this 
presumes the assumption of linear relationship between two-covariates# 
Therefore, positive correlation signifies a mutual increase between 
them, or a decrease of one as the other decreases# The simple correla¬ 
tion coefficient is defined as* 




J: (x,-x.) 




cc 


where and Xg are the co-variates. The arithmetic mean of each are 
represented by X-^ and respectively# The standard deviations of X^ 
and J -2 from the mean are ^x^ and <^X£ 8111(1 are the result of the formula* 
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where N is the number of co-variates under consideration# 

The line representing the linear regression of two variates is ex¬ 
pressed in the form# 


(X,-*.) ~ ~ (x 3 - X,) 

< 5 * 


a 
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"Regression of on Xg" merely states that an estimate of X^ may 
be obtained from a known value of X 2 through the use of the above formula 
in the form of X], a aX 2 + b. It should be emphasized that this equation, 
as so developed, is valid only for the estimation of X^ from X 2 and does 
not permit its use to estimate X 2 . 

The standard error of estimate, S, is defined as the mean of the 
differences between an observed X^ and the estimated values of Xj^ from 
above. It is a measure of the mean error involved when such an equation 
is used to estimate X^ from X 2 and, hence, may be considered as a standard 
for judging the reliability of the derived relation. S is, by definition, 
expressed as 




The calculation of multiple correlations involving three and four 
variates, as accomplished in this study, is also based on the assumption 
of linear regression. Such multiple correlation coefficients are ex¬ 
pressed by the relation: 


r x t -x 4 x 3 .... 
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where R is the determinantt 
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and Rj^ is the minor of rjjs 
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The standard error of estimate of X^ on the basis of Xg , X 3 , , 

Xy is defined as 
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And finally, the regression equation of on Xg, X 3 , — , Xy is 
given by t 



2* "Skill Score" computations^ ^ 

This is a statistical method for eliminating from a purely percent¬ 
age computation of the favorable distribution of a particular tetrachoric 
pattern, -that part of such a distribution which might result from chance* 
In effect, then, it is a percentage measure of the improvement of the 
given pattern over such a chance arrangement* 

The Skill Score, SS, is arrived at as follows: 

Given a tetrachoric distribution as: 

1 

1 /«. 

r. ! k 

I 

then 
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where E„ is the number of cases estimated to have been included in the 

v 

favorable quadrants by chance, and is defined as: 


cjx)' + 


It may easily be verified that, on the basis of the above develop¬ 
ment, a "pure chance" arrangement, i.e., Xj ■ Xg - Yj - Y 2 , would result 
in a Skill Score of SS « 0.0$* 
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